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Abstract 

Two  factors  that  influence  the  load  current  symmetry  of 
the  Atlas  pulsed  power  machine  are  analyzed  to  establish 
design  criteria.  Deviations  in  peak  current  among  the  24 
Marx  modules  providing  current  to  the  load  must  be 
limited  to  ±3.1%  to  achieve  the  goal  of  ±100  pm 
asymmetry  in  the  imploded  liner  radius.  Current 
asymmetry  may  also  be  caused  by  azimuthal  inductance 
variations  in  the  disk  transmission  line  leading  to  the  load. 
Mechanical  fabrication  tolerances  are  presented  that  must 
be  met  to  achieve  the  symmetry  goals. 


I.  INTRODUCTION 

The  Atlas  pulsed  power  machine  is  a  very  high  current, 
fast  discharge  capacitor  bank  presently  being  constructed 
[1]  at  Los  Alamos  National  Laboratory.  The  principal 
application  of  Atlas  will  be  to  implode  cylindrical  metal 
shells  (“liners”)  to  very  high  velocity  so  that  a  variety  of 
high  energy  density  hydrodynamic  experiments  [21  can  be 
performed.  Many  of  the  planned  experiments  require  that 
the  imploding  liner  accurately  maintain  its  circular  cross- 
section  during  the  implosion.  Deviations  from  circularity 
can  arise  from  azimuthal  variations  in  the  liner  density  or 
in  the  driving  current  density.  Experience  with  liner 
implosions  on  the  Pegasus  II  capacitor  bank  [3]  has 
clearly  demonstrated  that  material  and  fabrication 
tolerances  are  more  than  adequate  to  meet  the  symmetry 
requirements.  The  principal  concern  therefore  is  current 
symmetry. 

The  current  symmetry  goal  for  Atlas  derives  directly 
from  the  implosion  symmetry  goal.  Based  on  the  analysis 
presented  in  Section  II  we  expect  the  dominant  liner 
asymmetry  to  be  a  dipole  mode,  equivalent  to  a 
displacement  of  the  liner  axis  from  the  target  axis.  For 
this  asymmetry  mode  our  experimental  objectives  can  be 
met  if  the  deviation  from  an  ideal  circle  is  less  than  ±100 
pm  at  a  3:1  convergence  (roughly  15  mm  radius).  To 
achieve  this  in  a  typical  32  MA  implosion  to  12  km/s  the 
linear  current  density  at  the  entrance  to  the  load  region 
must  be  uniform  in  azimuth  to  ±0.35%. 

Current  asymmetry  in  Atlas  will  arise  from  two  primary 
causes,  variations  in  current  amplitude  among  the  24 
Marx  generators  supplying  current  and  deviations  from 
azimuthal  symmetry  in  the  disk  transmission  line  carrying 
current  to  the  load.  The  effect  of  current  variations  in  the 
Marx  generators  is  analyzed  in  Section  II  and  a  design 
criterion  is  derived  for  the  Marx  generator  performance. 
The  problem  of  transmission  line  asymmetry  is  addressed 


in  Section  III.  Guidelines  are  developed  for  the 
mechanical  tolerances  that  need  to  be  achieved  in  the 
fabrication  and  assembly  of  the  disk  transmission  line. 


II.  CURRENT  SOURCE  VARIATIONS 

Current  is  provided  in  Atlas  by  24  Marx  generators, 
each  delivering  1 .33  MA  with  a  risetime  to  peak  of  4.5ps. 
The  Marx  generators  are  housed  in  12  oil  tanks  located  in 
a  circular  pattern  approximately  7  meters  from  the  load. 
Current  is  carried  radially  inward  from  each  Marx 
generator  by  an  oil-insulated,  tri-plate  transmission  line. 
The  24  tri-plate  lines  join  to  a  common  conductor,  called 
the  “transition  section,”  at  a  radius  of  0.8  m.  From  0.8  m 
radius  to  the  load  at  0.045  m  radius  current  is  carried  by  a 
disk  transmission  line  that  has  both  flat  and  conical 
regions. 

Connections  to  the  transition  section  are  not  uniformly 
spaced.  The  pair-wise  mounting  of  the  Marx  generators 
and  the  need  to  provide  access  between  the  Marx 
generator  tanks  results  in  the  connection  geometry  shown 
in  Fig.  1.  The  connection  pattern  is  based  on  a 
fundamental  increment  of  9  degrees,  1/40  of  a  full  circle. 
The  two  transmission  lines  from  a  common  tank  are 
separated  by  9  degrees.  The  tanks  are  arranged  in  4 
groups  of  three.  An  extra  9  degree  space  separates  the 
tanks  in  a  group  and  there  is  a  wider  27  degree  separation 
between  groups.  Although  this  connection  pattern  is  not 
uniform,  it  has  a  high  degree  of  symmetry  including  a 


Figure  1.  Connection  geometry  of  24  Marx  generators  to 
central  transmission  line. 
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A.  Analysis  of  Radial  Current  Flow 

In  order  to  relate  variations  in  the  Marx  generator 
current  amplitude  to  the  symmetry  of  the  load  current  it  is 
necessary  to  account  for  the  smoothing  effect  of  the  radial 
current  flow  in  the  disk  transmission  line.  For  the 
purposes  of  this  analysis  the  disk  line  can  be  treated  as  a 
simple,  uniformly  spaced  pair  of  flat  disks  running  from 
0.80  to  0.045  m  radius.  This  is  a  reasonably  accurate 
approximation  for  a  conical  line  if  the  vertical  rise  is 
much  less  than  the  radial  dimension  as  it  is  for  Atlas. 

Two  further  approximations  can  be  made  that  expedite 
the  analysis.  First,  current  flow  can  be  calculated  in  the 
quasi-static  approximation.  The  time  scale  for  the  Atlas 
current  waveform  is  about  co''=3  ps.  The  transit  time  of 
an  EM  wave  through  the  disk  transmission  line  is  only 
Ar/0.65c  =  4  ns.  The  importance  of  displacement  current 
in  determining  current  flow  is  roughly  proportional  to  the 
ratio  of  these  time  scales,  —  1  O'3,  which  certainly  justifies 
the  use  of  the  quasi-static  approximation. 

The  second  approximation  concerns  the  relative 
importance  of  inductive  and  resistive  voltages  in 
determining  current  flow.  For  a  small  radial  section  of  the 
disk  line  the  ratio  of  inductive  voltage  to  resistive  voltage 
is 


^L  =  _i£_4,  (,) 

Vind  ^,<5(0  h  I 

where  p  is  the  conductor  resistivity,  8  is  the  skin  depth 
and  h  is  the  separation  between  conductors.  During  the 
current  rise  the  quantity  I/I-dot  ~  t  and  Eqn.  1  can  be 
further  simplified  to 


Eqns.  3  can  be  solved  by  separation  of  variables  to 
yield  the  following  expression  for  the  radial  current 
density. 


r  r 


n= 1 


C anrn  +-^-)cos(n<9)  + 
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( bnrn  +^-)sinOi0) 
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(4) 


The  azimuthal  current  density  has  a  similar  expression 
that  can  be  derived  from  Eqn.  4  using  Eqn.  3b. 

The  coefficients  {a,b,c,d}  are  determined  by  the 
boundary  conditions  at  the  inner  and  outer  radius  of  the 
disk,  r,  and  r0.  At  the  inner  radius  the  transmission  line  is 
short-circuited  by  the  load  requiring  that  Br  =  p^  Je  =  0. 
The  condition  Je(rj,t)  =  0  requires  the  following 
relationship  among  the  coefficients. 


b  r"  =  — 

n  1  rn 
ri 


(5) 


The  requirement  that  the  integral  of  rJr  over  all  angles  is 
equal  to  the  total  current  determines  ao  =  I,/2rc.  The 
remaining  coefficients  must  be  found  by  matching  the 
boundary  conditions  at  the  outer  radius.  Denoting  the 
current  input  from  the  \th  source  by  Ij(t)  and  the  angular 
position  where  it  is  connected  by  0j  the  radial  current 
density  at  rn  is  given  by 


1  24 

4(ro,0)  =  -£/,(O<5(0-0j)  (6) 

r°  M 

The  coefficients  \  and  bn  are  found  by  equating  Eqns.  4 
and  6,  multiplying  by  sin(m9)  and  cos(m0)  and 
integrating  over  0.  The  result  is 


2S0) 

nj 


At  early  time  the  skin  depth  is  small  and  inductance 
dominates  the  current  flow.  Toward  the  end  of  the  current 
rise  at  t=3  ps  the  skin  depth  for  an  aluminum  conductor  is 
0.3  mm.  For  a  disk  spacing  of  6.3  mm  the  voltage  ratio  is 
0.12  so  that  resistance  is  still  playing  a  minor  role  in 
determining  the  current  distribution.  For  the  remainder  of 
this  analysis  only  inductive  effects  will  be  considered. 

The  quasi-static  equations  for  the  magnetic  field  inside 
the  disk  transmission  line  are  VxB=V  B=0.  Using  the 
boundary  conditions  at  the  conductor  surface,  Bz=0,  Br  = 
p0  Jfl,  B0=  -p„  Jr,  these  equations  can  be  cast  in  terms  of 
the  conductor  surface  currents. 


drJa  dJ  r 

a)  - e- - L  =  0 

dr  30 


b) 


drJ  dj„ 
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24 

'Yjlj(t)  cos (mOj) 


7  =  1 
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(7) 


and  an  identical  expression  for  bm  with  sin()  replacing 
cos().  Finally,  the  desired  expression  for  the  current 
density  at  the  load  can  be  obtained  by  substituting  Eqns.  5 
and  7  into  Eqn.  4. 


JM,0) 


where 
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1  24  1  24 
cn  =—  £/,(f)cos(n0;)  Sn  =  —  2^Ij(t)s\n(ndj). 

1 1  7= i  1  7=1 

The  coefficients  (Cn,  Sn)  are  the  amplitudes  of  the  nth 
multipole  component  of  the  current  excited  at  the  outer 
radius.  Each  multipole  component  appears  at  the  load 
attenuated  by  the  ratio  (r,/r0)n.  Since  r/r„  =  0.056  for  Atlas 
it  is  apparent  that  only  the  lowest  order  components  will 
have  a  significant  effect  on  the  load. 


B.  Effect  of  Connection  Geometry 

Eqn.  8  can  be  used  to  evaluate  the  effect  of  the  non- 
uniform  connection  spacing  shown  in  Fig.  1  for  the  case 
where  all  of  the  input  currents  are  identical,  I,  =  I,/24.  For 
the  choice  of  axes  in  Fig.  1  the  0|  in  the  first  quadrant  are 
[1,5,7,13,15,19]  ji/40.  The  0j  in  the  other  quadrants  are 
equal  to  the  first  quadrant  angles  plus  7C/2, 7t  and  3rt/2. 

Table  1  displays  the  values  of  Cn  for  n=l  to  40.  For  the 
axes  chosen  all  Sn=0.  The  only  non-zero  coefficients  are 
those  where  n  is  an  integer  multiple  of  4  due  to  the  4-fold 
symmetry  of  the  connections. 


Table  1.  Current  distribution  coefficients  for  the 
case  where  all  input  currents  are  equal. 


n 

cn 

n 

cn 

4 

0.121 

24 

-0.167 

8 

-0.167 

28 

-0.513 

12 

0.513 

32 

0.167 

16 

0.167 

36 

-0.121 

20 

0 

40 

-1.000 

The  largest  asymmetry  induced  at  the  load  by  the  non- 
uniform  spacing  is 


other  and  so  forth  reduces  the  n=l  coefficient  to  a 
negligible  level. 

Excitation  of  the  n=2  and  higher  modes  is  not  a 
concern.  Consider,  for  example,  a  worst  case  where  2  of 
the  24  modules  located  opposite  each  other  fail  to 
discharge.  This  does  not  excite  n=l  but  it  does  produce  a 
C2-  0.08.  The  resulting  load  current  asymmetry  is  2x1 0‘\ 
well  within  the  design  requirement.  Since  n=2  and  higher 
are  so  strongly  attenuated,  only  the  n=l  mode  will  be 
addressed  below. 

Variations  in  current  amplitude  from  pulse  to  pulse  can 
be  represented  in  the  form 

24 

Jr(ro,d)  =  -f-]T(l  +  ej)S(0-0J),  (9) 

^ 1 O  ;  =  l 


where  the  Ej  represent  small  deviations  from  the  ideal 
current.  In  order  that  It  be  the  total  current,  the  Ej  must 
have  a  mean  value  of  zero.  Prior  to  operation  of  the  Atlas 
system  there  is  no  information  concerning  the  distribution 
function  for  the  £j.  Two  alternatives  will  be  considered.  In 
the  first  case  the  Ej  are  uniformly  distributed  over  an 
interval  ±£o.  The  second  case  is  a  normal  probability 
distribution  of  width  ±en.  Using  the  current  distribution  of 
Eqn.  9  and  restricting  the  analysis  to  n=l,  the  load  current 
distribution  is  given  by 


JM,9)  = 


In  r 


1  +  4  r'/  (Cj  cos  0  +  5,  sin  9) 


(10) 


where 


24 


24 


C\  ei  cos9J  £j  sin07 

7=1  7=1 


W  0.121  =  9.6  10~6, 


much  less  than  the  design  requirement  of  3.5x10 3. 


C.  Effect  of  Current  Amplitude  Variations 

The  24  current  sources  that  comprise  Atlas  will  not 
provide  identical  current  amplitudes.  One  may  expect 
both  reproducible  differences  due  to  manufacturing 
tolerances  in  the  components  and  fluctuating  differences 
due  to  small  variations  in  charging  voltage,  spark  gap 
impedance  and  so  forth.  In  general,  the  reproducible 
differences  are  not  of  much  concern.  Once  the  sources 
have  been  calibrated  it  is  relatively  simple  to  arrange  them 
around  the  disk  transmission  line  so  as  to  minimize 
excitation  of  the  n=l  mode.  For  example,  the  strategy  of 
placing  the  two  highest  current  sources  diametrically 
opposite  each  other  and  then  the  next  two  opposite  each 


The  peak  to  peak  variation  of  Jr  can  be  found  by 
locating  the  maximum  and  minimum  of  the  angular  part 
of  Eqn.  10.  The  final  result  is 


Eqn.  1 1  can  be  evaluated  statistically  to  place  bounds 
on  the  allowable  variation  in  current  amplitude.  Table  2 
presents  the  results  of  such  a  calculation  in  the  form  of  a 
probability  that  the  design  criterion  AJ/Jr  <  3.5x10-3  will 
be  exceeded  on  a  given  test  as  a  function  of  the  width  of 
the  distribution  function.  Each  probability  is  calculated 
by  selecting  a  24  component  £  vector  from  the  appropriate 
distribution  function,  calculating  C,  and  S,  according  to 
Eqn.  10  and  then  comparing  the  calculated  variation,  Eqn. 
11,  to  the  design  criterion.  This  calculation  is  carried  out 
N  times,  where  N  is  105  to  106.  The  total  number  of  times 
the  criterion  is  exceeded,  divided  by  N,  is  the  probability. 
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Table  2.  Probability  of  exceeding  the  load  current 
asymmetry  criterion  as  a  function  of  the  variation 
in  current  source  amplitude. _ 


Uniform 

Normal 

distribution 

distribution 

Eo 

Prob. 

£n 

Prob. 

.025 

.0005 

.0130 

.0005 

.028 

.003 

.0145 

.0017 

.031 

.009 

.0160 

.005 

.034 

.021 

.0175 

.011 

.037 

.040 

.0190 

.019 

.040 

.063 

.0205 

.033 

For  an  acceptable  failure  rate  of  1  in  100  the  amplitude 
variation  must  be  kept  to  less  than  ±3.1%  for  the  uniform 
distribution.  The  criterion  for  the  normal  distribution  is 
±1.7%.  This  appears  to  be  a  more  stringent  requirement. 
However,  the  normal  distribution  has  significant 
population  out  to  2a  (=  3.4%)  so  the  two  criteria  are 
roughly  equivalent. 

The  probability  of  exceeding  the  asymmetry  criterion  is 
a  sensitive  function  of  e.  A  20%  increase  in  e  causes  a  4- 
fold  increase  in  the  failure  rate.  For  this  reason  it  will  be 
important  during  commissioning  of  Atlas  to  measure  the 
amplitude  variations  and  identify  any  modules  with 
excessive  variations  so  that  the  source  can  be  identified 
and  eliminated. 


III.  TRANSMISSION  LINE  ASYMMETRY 

The  previous  section  demonstrated  the  beneficial  affect 
of  a  converging  disk  transmission  line  in  smoothing  out 
current  asymmetries  introduced  at  large  radius.  This 
section  will  address  the  complementary  effect,  the 
sensitivity  of  the  load  current  to  asymmetries  in  the  disk 
transmission  line. 

In  the  derivation  of  Eqn.  3  it  was  assumed  that  the 
separation  h  was  independent  of  9  and  r.  Allowing  h  to 
depend  on  6  introduces  an  asymmetry  in  the  disk 
transmission  line  inductance  that  can  perturb  the 
symmetry  of  the  load  current.  The  size  of  this  effect  can 
be  estimated  by  noting  that  for  small  perturbations, 
Ah/h«l,  the  magnetic  field  is  primarily  in  the  9 
direction  and  that  flux  conservation  requires  Be7i  = 
constant.  Combining  this  relation  with  the  boundary 
condition  for  B0  gives  for  Jr 

'-('’'“/'MO)  <12> 

Thus,  to  satisfy  the  requirement  AJ/J  <  3.5xlO'3  requires 
Ah/h  <  3.5xl0‘3.  This  is  a  very  stringent  mechanical 
tolerance.  For  a  typical  h  =  6.35  mm  (1/4”),  Ah  <  0.02 
mm  (0.001”). 

A  more  accurate  value  for  the  mechanical  tolerance 
could  be  obtained  by  repeating  the  earlier  analysis  and 
including  the  r  and  9  dependence  of  h.  However,  the 


solution  of  this  problem  is  quite  involved  except  for  a  few 
special  cases  and  was  deemed  unnecessary  for  the  present 
need.  Instead,  a  lumped  circuit  model  of  the  transmission 
line  was  formulated  and  solved  using  the  circuit  code 
SPICE.  The  lumped  circuit  model  is  shown  in  Fig.  2. 
Resistance  has  been  neglect  as  explained  above.  Fig.  2  is 
about  the  simplest  model  that  can  accurately  represent  a 
non-symmetric  current  flow.  The  angular  subdivision 
into  12  zones  allows  asymmetries  up  to  n=3  to  be 
investigated.  The  radial  zoning  is  based  on  a  logarithmic 
scale  so  that  each  zone  represents  about  the  same 
inductance.  The  choice  of  5  zones  allows  the  outer  zone 
to  represent  the  flat  portion  of  the  disk  line  while  the  inner 
conical  section  is  resolved  into  4  zones. 


Figure  2.  Lumped  element  model  of  disk  transmission  line 
used  to  calculated  effect  of  mechanial  variations. 

The  inner  edge  of  the  mesh  in  coupled  to  inductors  that 
represent  the  load  circuit.  Preliminary  calculations  using 
a  fixed  load  inductance  showed  that  the  current  symmetry 
was  sensitive  to  the  size  of  the  load  inductance.  For 
greater  accuracy  a  time  dependent  load  inductance  was 
added  to  the  model  to  approximate  the  behavior  of  an 
imploding  liner. 

Two  of  the  calculations  that  were  performed  with  this 
model  are  illustrative  of  the  conclusions.  In  Case  1  the 
axis  of  one  conductor  of  the  disk  line  is  displaced  an 
amount  e  with  respect  to  the  other.  This  produces  an  n=l 
asymmetry  in  the  conical  region  and  has  no  effect  on  the 
flat  region.  Some  typical  inductance  values  for  Case  1 
with  e=  0.07  h  are  shown  in  Fig.  2. 

The  calculation  was  performed  by  applying  identical 
current  pulses  of  the  form  I0  sin(cot)  to  each  of  the  outer 
nodes  and  measuring  the  current  flowing  in  each  of  the 
load  inductors.  The  asymmetry  for  Case  1 ,  defined  as  the 
maximum  difference  is  current  flow  at  peak  current 
divided  by  the  average  current,  was  3.4  x  10'2.  Thus,  in 
order  to  achieve  the  design  goal  of  AJ/J<3.5xl0'3  the 
offset  £<  0.007  h  or  Ah/h  <  0.01.  This  tolerance  is  about 
3  times  larger  than  the  estimate  derived  above.  The 
primary  reason  for  this  increase  is  the  smoothing  effect  of 
the  load  inductance. 
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For  Case  2  one  conductor  was  tilted  through  a  small 
angle  with  respect  to  the  other.  The  asymmetry  is  n=l 
again  but  the  principal  effect  is  at  large  radius. 
Summarizing  the  result,  the  allowable  deviation  at  the 
outside  radius  is  A hJh  <  0.04.  The  deviation  is  4  times 
greater  than  for  Case  1  because  Ah/h  is  small  near  the  load 
for  this  geometry  and  radial  convergence  can  smooth  out 
some  of  the  asymmetry  introduced  at  large  radius. 

Design  of  the  disk  transmission  line  is  proceeding  with 
the  objective  of  achieving  Ah/h  <  0.005  in  the  region  from 
4.5  to  15  cm  radius.  Achieving  the  highest  possible 
symmetry  in  this  region  is  doubly  beneficial.  It  avoids 
disturbing  the  current  symmetry  and  it  helps  to  smooth 
out  asymmetries  created  at  larger  radius.  Outside  of 
15  cm  radius  the  allowable  Ah/h  increases  linearly  up  to 
0.02  at  the  80  cm  outer  radius  of  the  disk. 
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